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Abstract. We consider the ensemble of n x n Wigner hermitian matrices 
H = {hik)i<e k<n that generalize the Gaussian unitary ensemble (GUE). The 
matrix elements h^i = hik are given by hik = {^ik + V~^yik)y where 

^ikyVek for 1 <: f- < k < n are i.i.d. random variables with mean zero and 
variance 1/2, yu = and xu have mean zero and variance 1. We assume the 
distribution of X£k , Vtk to have subexponential decay. In [3] , four of the authors 
recently established that the gap distribution and averaged fc-point correlation 
of these matrices were universal (and in particular, agreed with those for 
GUE) assuming additional regularity hypotheses on the xik^Vik- In 'Ji^ the 
other two authors, using a different method, established the same conclusion 
assuming instead some moment and support conditions on the x^k > Htk ■ In this 
short note we observe that the arguments of [3] and [7] can be combined to 
establish universality of the gap distribution and averaged fc-point correlations 
for all Wigner matrices (with subexponentially decaying entries), with no extra 
assumptions. 



1. Introduction 



This paper is concerned with the bulk eigenvahie statistics of large Hermitian 
Wigner random matrices, which we now pause to define. 

Definition 1 (Wigner matrices). Let n be a large number. A Wigner Hermitian 
matrix (of size n) is defined as a random Hermitian n x n matrix H — Hn with 
upper triangular complex entries hgk — n^^^^{xik + V—^Vik) for I < £ < fc < n 
and diagonal real entries hn = n^^^^xu for 1 < i < n, where 

• For I < £ < k < n, the xik, yik are i.i.d. copies of a real random variable x 
with mean zero and variance 1/2. 

• For 1 < £ < n, the xu are i.i.d. copies of a real random variable x with 
mean zero and variance 1. 

• x,x have subexponential decay, i.e., there are constants C, C" such that 
P(|a;| > t^) < exp(-t), P(|J| > t^) < exp(-t), for aU t > C . 
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We refer to x as the atom distribution of H. An important special case of a 
Wigner Hermitian matrix is the Gaussian unitary ensemble (GUE), in which x,x 
are Gaussian random variables with mean zero and variance 1/2, 1 respectively. 

Let Ai < . . . < A„ be the eigenvalues of a Wigner Hermitian matrix _ff„, and let 
Pn{xi, ■ . ■ ,Xn) be the associated (symmetric) density function, normalized to have 
total mass 1. For any 1 < k < n, the fc-point correlation function pi^^\xi, . . . ,Xk) 
is defined by the formula 

pli'^Xi, . . . ,Xk) / Pn{xi,. . . ,Xn) dXk+l ■ ■ - dXn- 



Here we follow the convention that the correlation functions are -normalized. 
Other conventions are also used in the literature, e.g. the correlation functions Rk 
used in [U 11] are related to plf^ via Rk = ^n"'k)\ P'^^ ■ 

The main results of this note are as follows. 

Theorem 2 (Universality of averaged correlation function). Fix e > and u such 
that —2 <u — e<u-\-£<2. Let k > 1 and let f : M'"' — > M &e a continuous, 
compactly supported function, and let H = Hn be a Wigner random matrix. Then 
the quantity 
(1) 

TT / / /("i'---'"fe)T r7vifcPnH^'+ — ^^r7T'---'"'+ — '^^^) dai . . . dakdu' 

converges as n —> oo to 

/(ai, . . .,ak) Aet{K {ai, aj))\,j^-^ dai . . . dak, 

where K{x, y) is the Dyson sine kernel 

, , sin(7r(x — v)) 

(2) K{x,y):^^ P 

and Psc{u) is the Wigner semicircle distribution 

(3) pUu):^^i^^u')T- 

Zn 

Theorem 3 (Universality of averaged eigenvalue gap). Fix u such that —2<u< 
2, choose k > 1 and s > 0, and let H be a Wigner random matrix. Then the 
expectation value of the empirical gap distribution 

An{u,s,e) = < j < n : Aj+i - A^- < ^— ; \Xj - u\ < e] 

satisfies 

r d"^ 

(4) lim lim EA„(u,s,e)= / — r det(l — /Cq) da, 

6->-0n->-oo Jq da 

where K-a is the integral operator on i^((0, a)) whose kernel is the Dyson sine kernel 
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We now briefly recall the previous partial results towards these theorems. In 
the case of GUE, the proofs of Theorems [21 [3] can be found in 5^. These proofs 
relied heavily on the explicit formula for the fc-point correlation functions in the 
GUE case. An important extension was made by Johansson [4], who established 
the above theorems in Gaussian divisible case when the matrix H took the form 
H + aV, where H is another Wigner matrix, is a GUE matrix independent of H 
and a > is an arbitrary fixed constant independent of n. This proof also requires 
the use of an explicit formula, however, its analysis is much more delicate than in 
the GUE case. By a simple rescaling, it is clear that Johansson's result extends to 
Wigner matrices of the form 

(5) = e-*/2ff + (l-e-*)i/2y, 

where H and V are as above, and the time t > is fixed independent of n. One 
can view H as the evolution of H under an Ornstein-Uhlenbeck process of time t. 

In [3] , four of the authors extended Johansson's result to obtain Theorems [2] and 
[3] for matrices of the form ([5]) with a short time t = n~^~^^ for any fixed 6 > 0. 
An important ingredient was the convergence of the density of the eigenvalues to 
the semicircle law ([3]) on 0{n~^'^^) scales which was previously established in [2]. 
Combining the result for Wigner matrices of the form ([5]) for time t — n^^^^ with a 
time-reversal argument, Theorems [2] and [3] were also obtained for Wigner matrices 
H with no Gaussian component, under the assumption that atom distribution 
of H has a probability density h{x) = e~f(^\ where g € C^(]R) grows at most 
polynomially fast at infinity and h{x) < Ce^"^!^'. For more details, we refer to [3l 
Section 1]. 

In the other two authors introduced a different way to study the local eigen- 
value statistics of random matrices, based on a variant of Lindeberg replacement 
strategy. This enables one to compare the statistics of a general Wigner matrix 
to those of a Johansson matrix, or any random matrix where the desired statistics 
is known, given that some moment condition is satisfied. As an application. The- 
orems m 13] were established in [7] under the additional hypothesis that the atom 
distribution x had vanishing third moment and was supported on at least three 
points. For more details, we refer to [71 Section 1]. 

In this note, we observe that one can combine the arguments from [3] with the 
arguments of [7] to eliminate the regularity hypotheses from the former and the 
moment and support hypotheses from the latter. The basic idea is simple and can 
be described as follows. The main result of [7] (see [71 Theorem 15]) allows one to 
compare the local statistics of two Wigner hermitian matrices whose atom variables 
f , C have the same third and fourth moments, namely E^'^ = E^''^ and E^** = E^'''. 
A closer look reveals that it is enough to require that the differences jE^^^ — E^''^| 
and |E^'' — E^'^'j are sufficiently small (in term of n). This allows us to compare the 
original matrix H with a Johansson-type matrix ^ with t ~ n^^^^ for some small 
(5, instead of comparing with a Johansson matrix with t = 0(1) as done in [7]. On 
the other hand, arguments from [31 guarantee that such a Johansson matrix has 
the desired statistics. 
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One defect of our method is that the statistics need to be averaged over a spatial 
scale e which is independent of n before universality is established. This appears to 
be a largely technical condition, arising from the lack of good concentration results 
for individual eigenvalues Ai for matrices of the form ([5]) (in particular, we do not 
know how to eliminate the vanishing third moment hypothesis from [71 Theorem 
29]). However, if we reimpose the hypothesis that the third moment of the atom 
distribution vanishes (i.e., Ex'^ = 0), we can remove the averaging in Theorem [21 
thus refining [71 Theorem 11] by removing the condition that the support must have 
at least 3 points (see Section [3]). 

2. Proof of theorems 

We now prove Theorem[21 We fix u, e, /, H; by a limiting argument we may take / 
to be smooth. We will always assume n to be sufficiently large depending on these 
parameters. Using the exponential decay of x and a standard truncation argument 
(see [3 Section 3.1]) we may assume that 

(6) X = 0(log°(i) n) 

almost surely. (Note that the Gaussian random variable N{0, 1/2) does not quite 
obey this bound, but can be modified on a set of extremely small probability to 
do so whenever necessary, so we may effectively assume that ([G]) holds in this case 
also.) 

We introduce the modified Wigner matrix 

H' :=e-*/2^+(l-e-*)i/V, 

where t := n^^^^ and is a GUE matrix independent of H. We will choose 
5 = 0.01 (say), but the proof below actually requires only S < 1/2. 

Observe that H' is another Wigner random matrix, with atom distribution x' 
given by 

:=e-*/2x + (l-e-*)i/2g 

where g = N[0, 1/2) is a Gaussian random variable independent of x. In particular 
we have the moment comparision bounds 

(7) I]{x'y = Ex^ 
for j — 0, 1, 2, and 

(8) E{x'y = Ex^ + ©(n-i+o-oi) 
for j = 3, 4. 

Let {p'nY'^'^ be the fc-point correlation function associated to H'. The first step is 
to show that Theorem [21 is valid for H' : 

Proposition 4. The quantity 

-L [ /(ai,...,a,)-4--(p;)WL'+^L-,...,u'+^f-)d«i...d«, 
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converges as n -> oo to 

/(«!,. . . , afe) det{K{at,aj))i ,^-^ dai . . . dak 



Proof. By the dominated convergence theorem, it suffices to show that 
(9) 

f{ai,...,ak)-, r\n:(Pn)^''\'^ ^ -Vt) - --dak 

f{ai, . . . ,ak)det{K{ai,aj))lj^^ dai . ..dak- 



lim 



for every fixed — 2 < u < 2. The dominated convergence theorem can be applied 
since the integral on the Ihs of (0) can be bounded by Ck E Nj , where Nj is the 
number of eigenvalues in an interval / of size C/n about u. The boundedness of 
ENj follows from [2] Theorem 5.1]. Originally this result was stated for atom dis- 
tribution with Gaussian decay, but the proof can be modified for random variables 
with subexponential decay and satisfying (|6]). 

Eq. ^ is proven in Proposition 3.1 of [3 . The proof is based on an explicit 
integral representation for the marginal {p[^)^^\ The only assumption used to prove 
([9]) in 3 is the asymptotic validity of the local semicircle law for the density of the 
eigenvalues of the Wigner matrix H (with no Gaussian part) on scales of order 
^-i+i5^ More precisely, if 

mn{z) — — trace , and msdz) ^ [ — ^^'^^ ^ 

n H — z J s — z 

denote the Stieltjes transforms of the empirical eigenvalue distribution and of the 
semicircle law, respectively, then we need to know [3j Eq. (3.9)] that for any k,5 > 
there exist constants c, C > such that 

(10) P( sup \mniE + iij)-msciE + iT])\>s]<Ce-''''^ 

\Ee{-2+K-2-K) J 

for n~-^^^ < < 1 and for all e > sufficiently small. 

Eq. (|10p has been established in [2] under the assumption of Gaussian decay of the 
atom distribution, that is assuming that Ee'^^ < cx), for some /3 > 0. To reduce this 
assumption to the desired sub-exponential decay assumption, we can use a result 
in [7] (see the paragraph following [T] Lemma 60]) which proved ([T0| under this 
weaker assumption with a weaker probability bound exp(— aj(log n)), where w(log n) 
denotes a quantity whose ratio with logn goes to infinity as n — >■ (X). One can also 
follow the ideas in Section 5 of [3], where it is shown that the Gaussian decay can 
be relaxed to exponential decay (at the expense of a slightly weaker probability 
bound). These ideas can be further refined to deal with sub-exponential decay. 

Following the proof of Proposition 3.1 in j3j, it is clear that the argument holds 
using this weaker bound on the probability as well. This concludes the proof of 
Proposition m □ 
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To conclude the proof of Theorem [21 it thus suffices to show that the expression 
^ only changes by o(l) when we replace H with H', where we use o(l) to denote 
an expression that goes to zero as n — cx). 

From the definition of pn , we can rewrite ([IJ as 

(11) - V EG(nA,,,...,nA,J 

n ^ — ' 

i<ii,...,ik<n distinct 

where 
(12) 

n /""^^ du' ( \ 

On the other hand, by 7, Proposition 61] (or Theorem 5.1]), taking advantage 
of the hypothesis ([6]), we see that, with probability 1 — 0(ri,~^°°) (say), 

#{1 < i < n : A, G /} < Cn\I\ 

for all intervals / C [u — e,u + e] with ]/] > \og'~' n/n, for some sufficiently large 
absolute constant C, C, where C can depend on u and e. Because of this, and the 
compact support of /, we can restrict the summation in (|TT|) to those ii, . . . ,ik for 
which ij = ii + 0(log'^*-^'' n) for all I < j < k, while only paying an acceptable 
error of o(l). In particular, there are now only 0(nlog*^^^'' rt) summands in pip . 
Similarly, using the convergence to the Wigner semicircular law (see e.g. [6]) we 
may restrict ii, . . . , ife to the interval [Sn, (1 — S)n\ for some fixed 6 > independent 
of n. 

A routine computation shows that G enjoys the derivative bounds 

lV^G(yi,...,yfc)l = 0(l) 

for all j/i, . . . ,2/fe and < j < 5, where the implied constant in the 0() notation 
can depend on u, e, /. 

Suppose first that that the random variables a;, x' matched to fourth order. One 
could then apply the four moment theorem in |7l Theorem 15] to conclude for each 
ii, . . . , ifc e [Sn, (1 — S)n], the quantity EG{nXi^ , . . . , nXii^) changes by 0{n^'^) for 
some absolute constant c > 0, uniformly in ii, . . . ,ife, when H is replaced by H' . 
Inserting this into (|11|) . we obtain the desired claim. 

In our situation, x and x' do not match exactly to fourth order, but only to second 
order. However, we have approximate matching to third and fourth order thanks 
to ([8|), and this is enough to recover the conclusions of the four moment theorem. 
Indeed, an inspection of the proof of [Zl Theorem 15] in |7i Section 3.3] reveals that 
the matching moments condition is used in only one place, namely to show that 
the expectation of an expression of the form 

F(0) + ^F{Q){h,k) + ... + ^V4F(0)(/itt, . . . , h,k) + 0(n-^/2+°(i)) 

only changes by at most nr'^^^-^ (say) when one replaces hu^ = •nT^/'^ {xgk + \/~ ^Vik) 
with h'^j. :— n^^/^{x£k+V-^yik), where I < i < k < n, and F : C — >■ E is a function 
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obeying the derivative bounds 



|V^F(0)| < Cn-^+°(i) 



for < J < 4. But the expectation of the first three terms of this expansion do 
not change when replacing hik with /i^^, thanks to (O, and the next two terms 



change by O(n-3+°(i)n3/2„-i+0-0i) ^nd C'(n-4+o(i)„4/2^-i+o.oi-) respectively, and 



the claim follows. 

Remark 5. The reader will observe that there is some room to spare here regarding 
powers of n. Indeed, a more careful examination of the argument reveals that one 
can sharpen Theorem [2] by not keeping e fixed with respect to n, but instead letting 
it decay a,s e — n~^^^^^ for any < S < 1/2. This should not be the optimal result, 
however; the correct scale should be 0(n^^+*) (and in fact one should not need 
to average in u at all). Unfortunately, our methods so far can only achieve this 
by putting more moment or regularity hypotheses on the distribution x, see for 
instance Section [S] 

We now sketch the proof of Theorem|31 The first step is to prove that the empirical 
gap distribution AJj of the matrix H' satisfies the equation Q. This follows from 
the same argument used in the proof of of Theorem 1.2 in [3]. The only difference 
is that the parameter e in ([4]) was replaced by (with some < 5 < 1) in [3]. 

Since the change of the density Psc(u') vanishes as e -> 0, the same proof applies to 
the current setting, using that the limiting function on the rhs of (jl]) is continuous 
in s. 

In the second step we deduce (|4]) for H by using that the expected values of observ- 
ables of the form (fTTj) are asymptotically the same for H and H' . By the exclusion- 
inclusion principle, A„ can be written as an alternating series of expressions in the 
form pi[) with functions G which are products of characteristic functions of inter- 
vals of length one or larger. These characteristic functions can be approximated by 
bounded smooth functions in L^-sense with arbitrary precision k. Controlling the 
alternating scries by its last term, letting first n — ^ oo then k — >■ and using the 
continuity in s of the limiting function, we obtain (|4j for H as well. 



If one assumes the additional third moment condition Ea:;'^ = on the moment 
condition, then it was shown in [Tj Theorem 29] that the eigenvalues with i g 
[Sn, (1 — d)n] obeyed the localization property 



with probability 1 — 0{n '-^) for any C, c > (and n sufficiently large depending 
on C, c, J), where t{a) was defined by the formula 

l-t{a) 



3. Concluding remarks 



A, = t{-) + 0{n 



l+c 




-2 
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Because of this localization, it is possible to eliminate the averaging over u' in 
Theorem [21 and to reduce the averaging in Theorem [3] to an interval of 0[n^) 
eigenvalues rather than 0(n), by combining the above the arguments with those 
used to prove [7j Theorems 9,11] ■ The net effect of this is to eliminate the hypothesis 
in those theorems that x is supported on at least three points. Thus, for instance, 
we can take x to be the Bernoulli variable that takes values ±1 with an equal 
probability 1/2. 

If the third moment condition from [71 Theorem 29] could be omitted, then one 
could similarly reduce the averaging required for Theorems [2l [3] even when the third 
moment was non-zero. To do so, however, would require either a relaxation of the 
moment conditions needed in ^ Theorem 15], or else an analogue of J, Theorem 
29] for the matrices H' studied here. It seems plausible that at least one of these 
two approaches could eventually be made to work, but this does not appear to be 
easily deducible from the existing literature. 

Another remark is that our theorems hold for a more general model of random 
matrices, where the real and imaginary parts of the entries are not necessarily 
independent. The reason is that 7], Theorem 15] (and its variant used here) do not 
require this assumption (see also the end of [71 Section 1]). 
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